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The configuration interaction (CI) method for calculating the exact eigenstates of a quantum- 
mechanical few-body system is problematic when applied to contact interactions between the parti- 
cles. In two and three dimensions, the approach fails due to the pathology of the Dirac (5-potential, 
making it impossible to reach convergence by gradually increasing the size of the Hilbert space. 
However, for practical applications this problem may be cured in a rather simple manner, by renor- 
malizing the strength of the contact potential, which must be diagonalized in a truncated Hilbert 
space. The procedure relies on the comparison of CI energies and wave functions with those obtained 
by the exact solution of the two-body Schrodinger equation for the regularized contact interaction. 
The rather simple scheme, while keeping the numerical procedures still elementary, nevertheless 
provides both cut-off-independent few-body physical observables and an estimate of the error of the 
CI calculation. 

PACS numbers: 67.85.-d, 31.15.ac, 34.50.Cx, 03.75. Ss 



I. INTRODUCTION 

Recent advances in the study of trapped ultra-cold 
atom gases [H 0, 01 have put great emphasis on the key 
role of two-body interactions in a whole class of phys- 
ical phenomena driven by particle-particle correlations 
Such systems represent an ideal laboratory for the 
investigation of many-body physics and the validation of 
many theoretical predictions. From the theory side, a 
major issue is the complexity of the atom-atom interac- 
tion potential, which can be hardly considered in full in 
theoretical simulations. On the other hand, very often 
only the low-energy scattering properties are relevant, 
which may be accounted for by simple model potentials. 
These model potentials may then be easily included e.g. 
in configuration interaction (CI) approaches to the few- 
body problem. The CI method displays unique features 
with respect to other numerical techniques, including an 
accurate treatment of correlation effects and access to 
both ground and excited states P|. 

A most simple and yet useful model potential is the 
contact potential @ , which describes the s-wave isotropic 
scattering by means of only one parameter, i.e., the scat- 
tering length a. However, in connection with quantum 
mechanics the use of pointlike interactions is mathemati- 
cally troublesome f^ , as also noticed in studies of bosonic 
few-body systems |^]. The reader may find a brief his- 
tory of contact potentials in Refs. H, [l0|, [O, to mention 
just a few. A brief review of the current applications of 
contact interactions in the field of ultra-cold atom gases 
is for example found in the introduction of Ref. IT^ . 

One major problem is that, in their simple form of 
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Dirac (5- functions, contact potentials do not provide any 
scattering in two and three dimensions [Tsl. fl3. [Tsl . 
This problem may be solved by introducing a proper 
regularization, which slightly enlarges the Hilbert space 
spanned by the possible eigenstates of the Hamiltonian 
This procedure, however, has subtle drawbacks, 
which are particularly relevant for the CI approach. The 
crux of the matter, as already noted by Huang 0, lies 
in the fact that the regularized operators may be non- 
Hermitian, which prohibits finding their eigenvalues by 
variational methods. In particular, the direct diagonal- 
ization of the Hamiltonian on a complete space provides 
only trivial (i.e., non interacting) solutions |6l. [14 1 . unless 
the class of allowed basis wave functions obey special 
(and often impractical) boundary conditions [l6l |. 

In this paper, we illustrate a rather simple alternative 
scheme to apply the CI method to the few-body prob- 
lem with contact interactions. In a truncated Hilbert 
space, the coupling constant of the interaction potential 
is renormalized according to the energy cut-off considered 
in the calculation. This relation is established for the ex- 
ample of two particles confined in a two-dimensional har- 
monic trap, for which the exact regularized solutions are 
available [1J| . We wish to point out that the CI diagonal- 
ization must be performed on a truncated Hilbert sub- 
space. In fact, the ground-state energy obtained from the 
CI solution of the two-body problem has no lower bound 
as the size of the basis employed in the diagonalization 
increases. Similar results were already reported by Esry 
and Greene ■ After truncation, the CI algorithm may 
correctly reproduce the exact energies and wave functions 
of both ground and excited states, except in a small range 
of the relative distance between the two particles, close 
to the origin. This agreement allows to map the coupling 
constant of the contact potential employed in the calcula- 
tion to the two-dimensional scattering length 020- This 
mapping, however, significantly depends on the energy 
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cut-off. 

In this paper, a relation between the couphng con- 
stant, the energy cut-off, and the scattering length is 
established, that straightforwardly may be used for CI 
calculations of few-body systems Jl8J. In spirit, this 
renormalization procedure is analogous to the ultravio- 
let renormalization adopted in similar contexts by means 
of diverse techniques, including space discretization [191, 
dimensional 11 1 and momentum-cutoff [ll|, [H, [20l. l2ll. 

regularization, as well as to other approaches 

The structure of the paper is the following: After ex- 
plicitely showing the failure of CI convergence in Sec. [Ill 
we discuss the exact solution of the two-body Schrodinger 
equation in Sec. IIIH which is based on the regularization 
of a pseudopotential operator. The key ideas are pre- 
sented in Sec. IIVI where we illustrate a successful strat- 
egy of comparing CI and exact results which provides a 
relation between the coupling constant of the contact in- 
teraction, the scattering length, and the energy cutoff of 
the calculation. After discussing the physical meaning of 
the truncation of the Hilbert space in Sec. |Vl we report 
the results of the renormalization procedure that are use- 
ful for actual CI calculations. Specifically, Sec. |VT] ilus- 
trates CI results for the two-body problem in the relative- 
motion frame, which hold for both fcrmions and bosons, 
while Sec. IVIII considers the case of Slater determinants, 
which is the most relevant for applications to N fermions. 
Conclusions are drawn in Sec. IVIIII 



II. THE TWO-BODY PROBLEM IN A 
TWO-DIMENSIONAL TRAP: FAILURE OF CI 
CONVERGENCE 

A simple, yet non-trivial example for investigating the 
nature of the contact interactions, is that of two inter- 
acting particles in a two-dimensional harmonic trap of 
characteristic frequency ujq. Throughout this paper we 
use as energy unit Tiluq and as length unit the oscillator 
length, io = {h/mujoy/^ . The two-body Hamiltonian, 
conveniently written in the center-of-mass and relative- 
motion coordinates R = {ri + r2)/2 and r = ri — r2, 
respectively, then reads as 
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where g is the dimensionless coupling constant in units 

The non-trivial part of H2, H,.^, pertains to the rela- 
tive motion, while the center-of-mass eigenstates of -ffcom 
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FIG. 1: Energy E of tlie two-body ground state in a two- 
dimensional harmonic trap vs number A'^i, of orbitals used in 
the CI expansion of the relative-motion wave function, plotted 
for different values of the coupling constant in the repulsive 
case, g > 0. Inset: same data for g = 0.1 with higher res- 
olution of the energy scale. The energy unit is the spacing 
hu!o between oscillator shells. The center-of-mass contribu- 
tion has been subtracted from the total energy, therefore in 
the absence of interaction E — 1. 



are simply Fock-Darwin orbitals in the R coordinate. 
The former may be solved by expanding the relative- 
motion wave function on a basis of Fock-Darwin or- 
bitals [13: 



*m(r) 



Nt-1 
n=0 



(r) 



(4) 



Here the 'Pnm('")'s are the eigenstates of Hrc\ when g = 0, 
and Nh is the number of orbitals employed in the di- 
agonalization. Moreover, if the relative-motion angular 
momentum is even (odd), the solution equally holds for 
spinless bosons and fermions with total spin zero (one). 
Here we investigate the lowest-energy singlet states of ei- 
ther two spinless bosons or two fermions of opposite spin 
1/2. The eigenvalue problem consists in the diagonaliza- 
tion of the matrix 



{n'\H,^i\n) ^ {2n+l)Sn^n' + 
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where \n) = (/?„(r) = ^Pnoif ) stands for the Fock-Darwin 
orbital in the r-space with n nodes in the radial direction, 
n = Q, . . . iVf, — 1, and to = 0, since the off-diagonal matrix 
element is zero for to 7^ 0. The explicit orbital form 
is ^„(r) = (27r)-i/2L„(rV2)e-'''/'^, where L„(z) is the 
Laguerre polynomial of order n. 

The ground state energy obtained from the numerical 
diagonalization of matrix (O is plotted in Figs. [1] and [5] 
as a function of the basis size iVb, for several values of 
5, both for the repulsive (Fig. [T|) and attractive (Fig. [J) 
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FIG. 2: Log-log plot of the energy E of the two-body ground 
state in a two-dimensional harmonic trap vs number Nt of 
orbitals used in the CI expansion of the relative-motion wave 
function, for different values of the coupling constant in the 
attractive case, g < 0. The center-of-mass contribution has 
been subtracted from the total energy, therefore in the ab- 
sence of interaction E = 1. 



case. From Fig. [T] (and the inset with higher resolution 
of the energy scale) it becomes apparent that the energy 
monotonously decreases as the basis size Ni, increases, 
not indicating any convergence for the basis sizes at hand. 
The failure of convergence is dramatically evident for g < 
-1 (note the log- log plot of Fig. 

This incapability of reaching convergence over a finite 
Hilbert subspace could be expected also from the fact 
that the off-diagonal Hamiltonian matrix elements are in- 
dependent from n [Eq. Considering the case of large 
interaction strength \g\, the Hamiltonian can be approx- 
imately written in matrix form as i?rei ~ -f + .9/(27r)l, 
where I is the identity and 1 is the Nh x Nf, matrix with 
all elements equal to one. For large values of Ni, the 
lowest eigenvalue is 1 and — \g\ / {2tt)Ni, for repulsive and 
attractive interaction, respectively, and it is clearly seen 
how the ground state energy diverges with N^, for at- 
tractive interaction [Tsj . This is also in agreement with 
the general expectation that the two dimensional Dirac 
(5-function does not provide scattering. Incidentally, we 
note that in the — well behaved — one-dimensional case 
the off-diagonal element is either zero or it decreases with 
n as - (-l)("+"')/2(W)-i/''. 



III. THE EXACT SOLUTION OF THE 
TWO-BODY PROBLEM 

An alternative approach to the solution of the two- 
body problem consists in introducing a properly regu- 
larized contact potential to replace the simple Dirac 6- 
function and then directly solving the Schrodinger equa- 
tion. This path was successfully followed by Busch and 
coworkers |14|, whose results we partly recall in this sec- 



tion referencing them as "exact". Later, Olshanii and 
Pricoupenko [28| provided a more general form of the 
pseudopotential, V^scudo, which should replace the sim- 
ple Dirac 5-function in i/rcb 
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where A is an arbitrary constant, a2D is the 
dimensional scattering length, A = e'''/2, and 
0.5772 ... is the Euler-Mascheroni constant. Here we ex- 
press a few remarks which are relevant to the following 
application of the CI method. 

The eigenstates of the regularized form of the relative- 
motion Hamiltonian, H^^^, with 
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are obtained by imposing that the wave functions, 'I'(t), 
written as a linear superposition of Fock-Darwin orbitals 
fnif) with unknown coefficients c„. 
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must solve the eigenvalue problem 

Hi:f^f{r)^E^{r). 



(8) 
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The condition ([9]) together with the requirement of nor- 
malization allow to analytically determine both the co- 
efficients c„ and the energy E, the latter through the 
equation 



^{l/2^E/2)^log{2/al^), 
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where ipi^) is the digamma function of argument z [2£ 
The wave function is 
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where v = {E — l)/2 and U{a,b,z) is the confluent hy- 
pergeometric function. For details of the derivation see 
Ref. M- 

Figure [3] shows the ground state energy E as a function 
of a2D, as obtained by solving Eq. (jTU]). We here disre- 
gard the lowest-energy solution of Eq. PH)) for < 020 < 
\/2, which is actually the true lowest-energy branch of the 
spectrum (not shown). The neglected branch refers to a 
deeply bound state whose energy indefinitely decreases 
as a2D 0^. The presence of such unusual bound 
state for repulsive interaction may be understood from 
the occurence of a pole in the upper plane of the com- 
plex wave vector of the scattering amplitude in free space, 
cf. Refs. and [s^. Throughout this paper we conven- 
tionally refer to the energy branch plotted in Fig. [3] as 
the lowest one. Such energy branch continuously evolves 
from the non-interacting ground state as a2D is var- 
ied, allowing for a one-to-one mapping between E and 
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FIG. 3: Exact energy E of the two-body ground state in 
a two-dimensional harmonic trap vs scattering length 020- 
The center-of-mass contribution has been subtracted from the 
total energy, therefore in the absence of interaction E = 1 
(dashed horizontal line). The dashed vertical line, for a2D ~ 
V2, separates the branch for repulsive interaction {E — > l""" 
as a2D ^0^) from that for attractive interaction {E ^ 1~ 
as a2D ^ +00). The energy unit is the spacing Tiuq between 
oscillator shells and the length unit is the oscillator length, 
£0 = {h/mcooy/^. 



a2D which will be extensively used throughout the pa- 
per. The energy E on the attractive side of the branch, 
for < a2D < +00, is always smaller than the non- 
interacting value E — 1, which is smoothly reached as 
a2D +00. For repulsion (0 < < \/2), E is larger 
than the non interacting contribution E ~ 1, approached 
for a2D 0+. 

In two dimensions, the interpretation of the scattering 
length is far from trivial: 020 is always positive and the 
cross section diverges in the zero-energy limit. Specif- 
ically, we may understand the result of Fig. [3] that E 
approaches the non-interacting limit as a2D — > +00, by 
recalling from scattering theory [SlJ] that a2D is the value 
of r for which the asymptotic form of the radial wave 
function r^/^^'(r) first becomes zero. In a short-range 
shallow attractive potential in three dimensions, the scat- 
tering length is negative for the lowest-energy scattering 
state, and it changes sign when the state becomes local- 
ized as the the depth of the potential increases. Exactly 
at the threshold depth, at which a bound state first ap- 
pears, the scattering length diverges. This explains the 
resonant behavior of Fig. [3] (a2D +00 as E — > 
since in two dimensions any arbitrarily weak attractive 
potential admits a bound state with radius exponentially 
large in the binding energy ^32]. As a2D is decreased on 
the attractive side, the binding energy becomes stronger 
{E is reduced) as well as the typical radius of the bound 
two-body complex decreases. In free space, there would 
be no lower bound for the energy [30|. The effect of 
the confinement potential is such that, as the scattering 
length is further decreased and it becomes comparable 



to the radius of the trap, 020 V^, the curve of the 
energy is folded into a different branch, which is actually 
the one we neglected. For further discussion see Refs. ISOl 
andil. 

The key point in the derivation of Eq. (|10p is that one 
is not allowed to interchange the limiting procedure con- 
tained in the definition ^ of Vpscudo, linrj.^o+ d/dr, with 
the infinite summation appearing in Eq. ([5]), X)n=o- 
fact, by doing so, the "regularizing" part of the pseu- 
dopotential, hmr^o+ ^og{AAr)d/dr, provides a null re- 
sult when applied to each addendum y„(r) of the sum, 
which is well behaved at the origin. In this case, the only 
possible solution of Eq. (fTO)) is the non-interacting one. 

For the same reason, the regularization of Ipsoudo is 
useless if the expansion ^ is finite, as it is of course the 
case for the CI algorithm. Thus, the regularization of the 
pseudopotential is irrelevant in this case. Equivalently, 
on any truncated Hilbert space, H^ci and H^°f of Eqs. ^ 
and ([7]) provide the same eigenvalues and eigenstates af- 
ter the identification 
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We conclude this section stressing that the physical 
observables of the two-body problem governed by H^^f do 
not depend on the numerical constant A appearing in the 
explicit form ([6]) of Vpscudo, as it may be straightforwardly 
checked e.g. in the derivation of Eq. (fTO|) . This fact is 
also illustrated by the solution of the scattering problem 
in free space associated to Vpscudo- From Appendix IXlone 
can see explicitely that the scattering amplitude does not 
depend on A as well as it coincides with the low-energy 
limit of any short-range potential. 

IV. COMPARISON BETWEEN EXACT AND 
CI RESULTS 

Let us now investigate whether it is possible to per- 
form a CI diagonalization in a truncated Hilbert space, 
corresponding to a physically sensible cut-off, in order to 
obtain results which are nevertheless cutoff-independent. 
To accomplish this we suggest to employ the degree of 
freedom associated to the choice of A in the expression 
p2p for g in order to renormalize the coupling constant 
itself. The physical requirement validating the procedure 
is that one obtains the same CI energies and wave func- 
tions for different values of the energy cut-off, as well 
as that the CI data compare well with the exact results 
available for the two-body problem. 

A. The two-body ground state 

CI energies are obtained by performing the diagonal- 
ization by fixing the value of the coupling constant g. 
The result depends on the number Nf, of orbitals used 
as a single-particle basis, as shown in Fig. [D By using 
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FIG. 4: (Colour) CI energy E of the two-body ground state 
in a two-dimensional harmonic trap vs coupling constant g, 
for different numbers A^;, of orbitals used in the CI expansion 
of the relative-motion wave function. The center-of-mass con- 
tribution has been subtracted from the total energy, therefore 
for g — — absence of interaction — E = 1. 

just one orbital, Nb = 1, one recovers the result of per- 
turbation theory at first order in g, therefore E depends 
linearly on g (cyan straight line in Fig. [4|) . By increas- 
ing Nh, E is reduced in a non-linear fashion, as it could 
be expected for a generic variational calculation (Nb in 
Fig. |4] increases by going from the green to the red curve). 
Since in the present case of contact interaction the diago- 
nalization does not converge as is increased, one needs 
to show that the same CI observables may be obtained 
for different values of N},. 

The first step in this empirical proof of principle is 
illustrated in Fig. [51 displaying the root mean square ra- 
dius for the relative motion, ri-.m.s., as a function of the 
energy E of the two-body ground state. This is the sim- 
plest quantity providing information on the ground state 
wave function, 

{poc -\ 1/2 J- ^oo 'I -1/2 

dr\^ir)\\H dr\^ir)\\j 

Obviously, repulsive (attractive) interaction implies an 
increase (a decrease) of both energy and radius with re- 
spect to the non-interacting values, E = 1 and rr.m.s. = 
V2, respectively. Figure [5] shows both exact and CI 
results: the possible ground states, represented in the 
(i?, rr.m.s.) plane, form a geometrical locus parametrized 
by the allowed values of either the scattering length 
for the exact case [black curve, cf. Fig.[3]and Eq. (fTO|) ]. or 
the coupling constant g for the CI data (coloured curves, 
cf. Fig.|4]). The black curve of Fig. [5] represents the exact 
data for the energy branch of the spectrum reported in 
Fig.[3l whose lower and upper bounds are E = —1.923264 
and E — 2.008166, respectively. The coloured curves are 
obtained from CI diagonalization for a broad range of 
the couphng constant, — 5 < .g < 30. In Fig. O the 
non-interacting case is represented by the point of co- 
ordinates {E, rr.m.s.) = (l,\/2), identified by the inter- 
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FIG. 5: (Colour) Root mean square radius, rr . m.s., VS energy 
of the two-body ground state, E. Both exact (continuous 
black curve) and CI data (coloured curves) are displayed: the 
former are obtained by varying the two-dimensional scatter- 
ing length, < a2D < oo (cf. Fig. [3]), the latter by varying 
the coupling constant, —5 < g < 30, for a given number of 
orbitals. A*";,, used in the CI expansion of the relative-motion 
wave function (cf. Fig. [4]). The black vertical dashed line is a 
guide for the eye: its intersection with the continuous black 
curve identifies the non-interacting ground state. CI data on 
the left (right) hand side correspond to attractive (repulsive) 
contact interaction. The center-of-mass contribution has been 
subtracted from the total energy. In the absence of interac- 
tion, the energy is _E = 1 and the root mean square radius is 

^r.m.s. — 



section of the vertical dashed line with the continuous 
black curve. The left (right) neighborhood of the non- 
interacting point corresponds to either the coupling con- 
stant g ^ 0^ (g 0+) in the CI case (cf. Fig. [4]), or the 
scattering length a2D ^ +oo (a2D ^ 0+) in the exact 
case (cf. Fig. [3]). 

From Fig. [5] we see that all CI curves, shown for dif- 
ferent values of Nf, , remarkably overlap in some region of 
the (E, rr.m.s.) space, centered around the non- interacting 
point (-E,rr.m.s.) — (1, V2)- The actual range of the E- 
axis along which CI and exact data overlap depends on 
the values of both Nj, and g: it is just the point E = 1 
for Nf, = 1 (cyan line), corresponding to the result of 
first-order perturbation theory that the wave function 
is not affected by the perturbation, and then it mono- 
toneously increases with Nb up to a maximum interval 
approx. 0.5 < E < 1.5 for Nb = 10 (red curve). 



B. Two-body excited states 

By comparing to the exact two-body solution p^ . 
Fig.[5]in principle allows to estimate the relation between 
g and a2D once E and Nb are fixed, as well as to assess 
the error associated to the CI diagonalization. However, 
there is still some arbitrariness in the comparison be- 
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TABLE I: Comparison between exact and CI energies for 
both ground and excited states, for attractive (2nd and 3rd 
column) and repulsive (4th and 5th column) interaction, re- 
spectively. Exact and CI data are linked by matching the 
ground state energies, providing (020,5) = (11.71,-1.855) 
and (0.3294,30.78) for E = 0.5, 1.5, respectively, where 020 
is obtained by solving Eq. (I10|l and g is obtained by diag- 
onalizing Eq. ([5]) with A^i, = 13. Excited state energies are 
then obtained by either solving Eq. (|10p with a2u fixed or by 
diagonalizing Eq. ((U with g fixed. 



tween CI and exact data, since we could have chosen as 
weh observables other than rj-.m.s. in order to character- 
ize \I'(r). Therefore, as a second step of our proof of 
principle, in the following we straightforwardly compare 
both energies and wave functions of a few excited states. 
Such criterion is obviously much more stringent than the 
comparison of Sec. IIV Al 

As an example, we fix the CI basis size {Nt, = 13) and 
consider two cases, corresponding to attractive and repul- 
sive interaction, respectively. The corresponding values 
of the coupling constant g are obtained by matching the 
CI ground state energies with the chosen reference val- 
ues. For E — 0.5 and E = 1.5 we obtain g — —1.855 
and 30.78, respectively. By fixing g we compute CI ener- 
gies and wave functions for the excited states, which arc 
compared to the analogous exact quantities derived for 
fixed scattering length a2u through Eqs. (fTU|) and pT|) . 
respectively. The scattering length is deduced by solving 
Eq. Uni) for E = 0.5 and 1.5, giving a2D = 11-71 and 
0.3294, respectively. 

Table [J compares energy levels up to the third excited 
state, showing that the worst relative error is 7 parts 
per thousand (for the 3rd excited level in the repulsive 
case) . This very good agreement is confirmed by the wave 
function analysis, reported in Figs. Elja) and (b) for at- 
tractive and repulsive interaction, respectively. The plots 
show the exact (red curves) and CI (black curves) radial 
probability densities of finding two particles at distance 
r, up to the third excited state. Such quantity, defined 
as r|^'(r)| and normalized to one, allows to quantify 
the discrepancy between exact and CI predictions. We 
see in both attractive [Fig.IB^a)] and repulsive [Fig.[6jb)] 
cases that the overlap of the probability densities as well 
as the agreement regarding node locations is almost per- 
fect far from the origin, becoming progressively worse as 
r — > 0. Such deviations, however, remain small in the 
whole range of r, so the overall matching between exact 
and CI results is excellent, consistently with the data of 
Table |T1 A closer inspection of the wave function behav- 
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FIG. 6: (Colour) Radial probability density r |^'(r)| of find- 
ing two particles at distance r, from both exact (red lines) 
and CI (black lines) calculations, for attractive (a) and re- 
pulsive (b) interaction, respectively. The panels display the 
ground (GS), first- (1st), second- (2nd), and third-excited 
(3rd) states obtained from both exact and CI calculations. 
Here '^(r) is the relative-motion wave function normalized 
such that J drr\^{r)^ = 1. For the exact calculations the 
values of 020 = 11.71 and 0.3294 were used for attractive (a) 
and repulsive (b) interaction, respectively, while the coupling 
constant used in the CI diagonalization was g = —1.855 and 
30.78, respectively, for a single-particle basis of A'^;, = 13 or- 
bitals. The corresponding energies, up to the third excited 
state, are displayed in the 2nd and 3rd (4th and 5th) columns 
of TableUfor attractive and repulsive interaction, respectively. 



ior close to the origin, as discussed below, reveals signifi- 
cant qualitative differences between exact and CI results, 
which however do not affect the substantial agreement 
between wave functions in the whole space. 

To summarize, we have shown so far that the CI 
method predicts energies and wave functions that com- 
pare well with the regularized ones, provided one works 
in a certain region of the {E^r-c ^^ s.) plane. This paves 
the way to the renormalization of the coupling constant 
g and to its actual usage in CI calculations. In the re- 
maining part of the paper we will consider in detail such 
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1 1.2 1.4 1.6 1.8 

E 



FIG. 7: (Colour) As Fig. [5] above, but zooming in on the 
region 1 < E < 2. The separation between two consecutive 
circles in each CI curve corresponds to an increment of the 
coupling constant Ag — 5. 





S = 0.5 (3<0) 


E = 1.5 (g > 0) 


B = 1 (<? = 0) 


2 


1.58 


3.54 


1.41 


3 


1.21 


2.57 


1.08 


4 


1.00 


2.14 


0.912 


5 


0.86 


1.86 


0.803 


6 


0.77 


1.67 


0.726 


10 


0.53 


1.21 


0.552 


13 


0.41 


0.99 


0.481 



TABLE II; Real-space cut-off, Vc, for a given number of Fock- 
Darwin orbitals, A^;,, employed in the CI diagonalization of the 
two-body ground state wave function 'I'ci('~) in the relative- 
motion frame of reference (see text). Data of second and 
third columns correspond to energies E — 0.5 and E — 1.5, 
being the largest roots of the equation D{rc) = -Dthrcshoid, 
with -Dthreshoid = 0.02 and 0.01, respectively. The last column 
reports the first zero of the A^th non-interacting orbital along 
the radial direction as a reference. 



perspective. 



V. ROLE OF CI TRUNCATION OF THE 
HILBERT SPACE 

In this section we discuss the physical significance of 
the CI truncation of the Hilbert space. Going back to 
Figs, m and EJ we note that, once a positive (negative) 
value of g is fixed, the corresponding energy obtained 
through the CI method decreases when Ni, increases. 
This may be seen in Fig. [5] e.g. from the displacement 
of the hmits of the different CI curves for increasing val- 
ues of Ni,. Such hmits correspond to g — 30 {g — —5) 
at the right (left) end of the curves. A magnification of 
Fig. [5] for a larger set of data in the region g > shows 
that the scaling of E with g is non-trivial. In Fig. [7] the 
separation between consecutive circles in each CI curve 
corresponds to an increment in the g-space of Ag ~ 5. 
As g increases, the circles tend to cluster around a point 
of accumulation, i.e., the CI energy E tends to an upper 
bound as g "> oo, the larger Nb, the smaller the limiting 
value. For g < 0, the situation is complementary (data 
not shown): s-— ooasg-^— oo, the larger N^, the 
larger \E\. By extrapolating these results to the limit of 
completeness of the Hilbert space, A^b ^ oo, it is clear 
that E' — > 1 for any fixed value of g > and E ^ — oo for 
g < 0- This behavior supports the fact that the contact 
interaction does not provide any scattering in the free 
two-dimensional space. 

The above analysis for repulsive interactions implies 
that a truncated, smaller Hilbert subspace allows to ac- 
cess a broader range of the a2D-space, which may be de- 
rived through Eq. ifTO]) or Fig. [3] linking E and a2D- How- 
ever, this criterion must be balanced by the accuracy of 
CI energies and wave functions, which decreases with the 
size of the Hilbert space, as we now discuss. 



In order to assess the accuracy of CI calculations, we 
define the real-space cutoff Tc as the lower bound of the 
interval (re, +oo) within which CI and exact wave func- 
tions overlap. Specifically, define the integral function 
D{r) as the deviation of CI radial probability density 

r|^'ci(f)| from the exact quantity r j^* exact (j^)! , inte- 
grated from r to oo: 



D{r) 



dr' 



l*Ci(r')r 



t(r')r 



(13) 



Since the wave functions are separately normalized to 
one, < D{r) < 1. Clearly, D{r) ^ as r ^ -|-oo since 
the wave functions vanish at infinity. As r is reduced, 
D{r) remains zero as long as ^cii'r) and \E'exact('') over- 
lap. As r is further decreased, ^ciir) starts to deviate 
from ^'exact(''), so D{r) reaches a global maximum at 
some finite value of r and then decreases again because 
D{r) for r ^ 0. We define Tj. as the position at 
which D{r) first assumes the value of r>threshoid as r is 
reduced from infinity. 



threshold- 



(14) 



For given Ni,, the value of g employed in the CI di- 
agonalization is chosen such that the CI energy matches 
the reference value of the exact energy. This approach 
is the same as that used for Table U and Figs. [DJa) and 
(b). The values of Tc as a function of for the ground 
states corresponding to the top left panels of Figs. [H^a) 
(-Dthrcshoid = 0.02) andEfb) (Ahrcshoid = 0.01), respec- 
tively, are collected in Table HIl The value of Tc manifestly 
decreases as Nf, increases, the larger \g\ the larger its mag- 
nitude. As a reference, we also report the first zero of the 
Afcth Fock-Darwin orbital along the radial direction, ob- 
tained as the smallest root of the equation (pN^^i{r) = 0. 
All the data show similar trends and may be fitted by 
a power law of the type cx n!^ , where P « —0.7 for 
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FIG. 8; Square modulus of the relative-motion wave function 
of the two-body ground state, |'I'(r)|^, vs radius, r, and as a 
function of the number of orbitals, Ni,, used in the CI diag- 
onalization. The value of the coupling constant, g < 0, for 
each subspace fixed by A'';, has been chosen according to the 
renormalization procedure described in the main text. One 
finds g = -2.40, -2.10, -1.85, for Nt = 3, 6, 13, respectively. 
The continuous line is the reference square modulus of the 
exact wave function of the top left panel of Fig. [6ja) , whose 
energy is E — 0.5. The wave function normalization is such 
that Jdr l*(r)l^r = 1. 

E = 0.5 as well as for E = 1.5, and (3 « —0.6 for the 
non-interacting case. Therefore, since the value of iVf, is 
the cut-ofF of the Hilbert space on the energy axis, Ec, 
and roughly N}, cx 1/r^, we may conclude that the cut- 
ofF values in real and energy spaces are complementary, 
Ec oc Here we assume that the cut-off values re- 

spectively in momentum space, Pc, and in real space, Tc, 
are related by pc cx l/vc- 

It is also interesting to investigate in more detail the 
deviation of the CI wave function from the exact one 
close to the origin. To this aim we r^lace the plots 
of the radial density probability r |5'(r)| of the top left 
panels of Figs, ^sl) and [6jb), respectively, with those 
of the square modulus |^'('')| , shown as a function of 
r in Figs. [8] and [H respectively. While the former is 
normalized to one and goes to zero as r — > 0, the latter 
may display possible singularities at the origin. 

Figures M {g < 0) and [9] (g > 0) compare the exact 
wave function (continuous line) with those obtained by 
diagonalization for various basis sizes Nf,. The square 
moduli of the exact wave functions present a logarithmic 
singularity at the origin, originating from the behavior of 
the hypergeometric function of Eq. pT|) : 

*exact(r-) ~ jrp^ [logr-^ + Tpi-i')] for r w 0, 

where A is a normalization constant. Note the appear- 
ance of a node in the exact wave function of Fig. [9l due 
to the occurrence of a deeply bound state at much lower 
energy (cf. discussion in Sec. IIIip . Any CI expansion of 




FIG. 9: Square modulus of the relative-motion wave function 
of the two-body ground state, |'I'(r)|^, vs radius, r, and as a 
function of the number of orbitals, A^;,, used in the CI diag- 
onalization. The value of the coupling constant, g > 0, for 
each subspace fixed by Ni, has been chosen according to the 
renormalization procedure described in the main text. One 
finds g = 6.00, 10.0, 30.8, for Nt = 3, 6, 13, respectively. The 
continuous line is the reference square modulus of the exact 
wave function of the top left panel of Fig. [6jb) , whose energy 
is i5 = 1.5. The wave function normalization is such that 
/dr|*(r)|V = 1. 

exact (f) over a finite set of basis functions regular at the 
origin is unable to reproduce the logarithmic singularity, 
therefore any sensible comparison between CI and exact 
results must regard only the form of ^oxact (r) sufficiently 
far from the origin. On the other hand, as Nh is increased 
(going from dotted to dashed-dotted and dashed curves 
in Figs. [5] and [5]), the overlap between exact and CI wave 
functions manifestly increases, as well as Tc diminishes, 
as shown in Table |TT1 The presence of a singularity for 
^cxact(f) is not important for practical calculations as 
long as its contribution to the radial density probability 
r |^cxact('')|^ is small, as demonstrated in Figs.lHJa) and 
Hb). 

The above scenario is implicit in the usage of the short- 
range pseudopotential as well, since that by definition is 
introduced in order to reproduce the asymptotic behav- 
ior of the scattering states (cf. Appendix On the 
other hand, similar concerns regard interaction poten- 
tials other than contact interactions. For example, CI 
diagonalizations employing simple orbital bases are un- 
able to reproduce cusps in the wave functions originating 
from Coulomb interactions 1331 1. 



VI. RENORMALIZATION OF THE COUPLING 
CONSTANT 

Let us now numerically investigate the functional rela- 
tion between, on one side, the value of the coupling con- 
stant g and the size of the Hilbert space Nf, actually em- 
ployed in the CI diagonalization, and, on the other side. 
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the physical observable describing the short-range scat- 
tering, i.e., the two-dimensional scattering length a2Dj 
whuch is linked to the exact energy E via Fig. [3l This 
mapping may be seen as a renormalization of the cou- 
pling constant g through the free parameter A appearing 
in the expression (|12|) , since the value of g corresponding 
to a certain value of a2D depends on Ni, as well. 

The rather simple idea, based on the plots shown in 
Figs. [Hand [7l is now the following: The value of g to be 
linked to a certain value of the exact energy — say E — 
is found by the intersection of the vertical line E = E 
with the CI curve corresponding to a given size Nh of 
the Hilbert space, since each point of the CI curve is 
associated to a specific value of g. The vertical distance 
between the CI and the exact black curve provides an 
estimate of the error of the CI calculation. Note that for 
two particles the value of E may be easily mapped to a2D 
through Eq. (fTU)) . 

Figures fTOTa) and (b) for repulsive and attractive inter- 
actions, respectively, show the key results of the paper. 
The plots allow for identifying either the value of the 
two-body exact energy, E (left vertical axis), or, equiva- 
lently, the two-dimensional scattering length, 020 (right 
vertical axis), as a function of the values of the coupling 
constant, g (horizontal axis), and of the size of the Hilbert 
space, Nb [coloured curves, going from Nt ^ 2 (orange) 
up to Nb = 10 (black)], employed in the CI diagonaliza- 
tion. The error associated to the comparison between 
CI and exact results is provided by the vertical error 
bars of the curves of Figs. fTOTa) and (b). Such error 
is estimated as follows: From plots like Figs. [5] or [7] we 
associate to a given value of the energy, E, a specific 
value of Nb and g: The latter is obtained by the inter- 
section of the vertical line E ^ E with the CI curve 
for Nb- The vertical distance between the CI and the 
exact curve provides the error for rj-m s., Arr.m.s.- Al- 
ternatively, from the intersection of the horizontal line 
''r.m.s. = ^r.m.s. with the CI and exact curves, one could 
obtain an estimate for the error on the energy, AE. If the 
exact and CI curves are sufhciently close, we may assume 
that Arr.m.s./Ai? « 9rr.m.s./9i?, where the latter is the 
slope of the exact curve. Therefore, from the knowledge 
of both Arr.in.s. and 9ri-.in.s./5i?, one obtains the estimate 
for AE that is actually plotted in Figs. [TUT a) and (b). 

The main features of the plots are consistent with the 
findings of Sec. |Vl First, the accuracy of the calculation 
is maximum close to the non-interacting point, g — s- 0, 
and it increases with Nb (the error bars are smaller if 
Nb is larger, once E is fixed). Second, if Nb increases, 
the energy depends more strongly on the strength of the 
contact interaction, g, for attractive interactions, while 
the opposite holds for repulsive interactions. 

We would finally like to remark that the procedure 
reported in this section is similar in concept to the re- 
cent renormalized mean-field theory developed by von 
Stecher and Greene, where the scattering length used in 
a Hartee-Fock calculation is renormalized by mapping 
the self-consistent energy of two particles to the exact 





FIG. 10: (Colour) CI energy of the two-body ground state, 
E, vs. coupling constant, (/ > (a) and g < (b), for dif- 
ferent numbers of orbitals, Nt, used in the CI expansion of 
the relative-motion wave function. The vertical error bars are 
estimated according to the procedure described in the main 
text. The right vertical axis shows the two-dimensional scat- 
tering length, a2D, obtained through Eq. (|10p . The center-of- 
mass contribution has been subtracted from the total energy, 
therefore the non- interacting energy is i? = 1. 



one 



VII. SLATER-DETERMINANT VS 
RELATIVE-MOTION BASIS 



So far we have considered only two interacting par- 
ticles, and for this special case we have uncoupled the 
relative motion from the center-of-mass oscillation. One 
of the advantages is that the results for two spin- 1/2 
fermions forming a singlet are the same as for two spin- 
less bosons. However, for iV > 2 this transformation 
turns out to be cumbersome for any efficient implemen- 
tation of the CI algorithm. In fact, in typical codes, 
the A^-body Hilbert space is spanned over a basis of ei- 
ther Slater determinants for fermions or permanents for 
bosons. It is therefore essential to verify that the results 



10 





1 


(a) / 




G— G 2 shells 


1.7 


o-o N^ = 2 




o—o 4 shells 






G— G N^ = A 






o—o 10 shells 




E 1 R 


c^A/^=10 










1.5 


1 





1 1.1 1.2 1.3 1.4 

£ 



1.4 


(b) ' ^ 


"?1.2 




E 


^^.--e-'^O;^ G— 2 shells 




_______e ''O::^^ C3-oA/^ = 2 




^Pji^-'^^^^^ ^=>—^ 4 shells 




^^^^^5^^^^ G— o A/^ = 4 


1 


3=^*="*^^^^^ G— o 10 shells" 




(^A/^=10 

1 



0.25 0.5 0.75 

£ 



FIG. 11; (Colour) Root mean square radius, rr.m.s., vs en- 
ergy of the two-fermion CI ground state, E, computed by 
varying g > (a) and g < Q (b) in different Hilbert sub- 
spaces employed in the CI diagonalization. Green, cyan, and 
red curves are obtained from a Slater-determinant basis built 
including Fock-Darwin levels up to 2, 4, 10 energy shells, re- 
spectively. Blue, magenta, and black curves are obtained in 
the relative-motion frame, including up to Ni, = 2, 4, 10 or- 
bitals, respectively. The separation between two consecutive 
circles in each CI curve corresponds to an increment of the 
coupling constant Ag — 1, with (; ^ 0^ as i5 1"*" (a) and 
p ^ 0~ as -B ^ 1~ (b). The center-of-mass contribution 
has been subtracted from the total energy. In the absence of 
interaction, E = 1 and rr.m.s. = V^. 

reported in the previous sections still hold for these new 
bases. 

To this aim we focus again on the case of two fermions. 
The Hamiltonian H2 of Eq. ([1]) , now written in terms of 
the standard coordinates ri and r2, reads 

^2 = E {-l^' + ^^^^ 

In order to apply the renormalization scheme to a CI cal- 
culation performed on a basis of Slater determinants, we 
start again from the relation rr.m.s. vs E (cf. Figs. [5] and 



[T]). To evaluate the square modulus of the relative- motion 
wave function, |\I'(r)|^, needed to compute the root mean 
square relative distance, rr.m.s. , we first compute the pair 
correlation function P(ri,r2): 

1 ^ 

(16) 

where the average, (...), refers to the CI ground state, 
and N = 2. By integration of P{ri,r2) over the center- 
of-mass coordinate, R, one is left with the probability 
density in the relative coordinate, |^'(r)|^: 

\^{r)f ^2TrJdRPir/2 + R,-r/2 + R), (17) 

where \I'(r) is circularly symmetric. 

Figures [TlT a) and (b) compare the CI curves, rr.m.s. 
vs E, for the two possible choices of the CI bases, ei- 
ther the Slater determinants (green, cyan, red curves) or 
the relative-motion Fock-Darwin orbitals (blue, magenta, 
black curves), in both cases of repulsive [g > 0, Fig. fTTT a)] 
and attractive [g < 0, Fig. fTTT b)] interactions, and for dif- 
ferent basis sizes. Since a major purpose of the present 
scheme is the application to systems with TV > 2, we 
limit ourselves to small basis sizes (maximum ten shells 
of Fock-Darwin orbitals). The CI code DONRODRIGO 
,35i] was used for the calculation on the basis of Slater de- 
terminants. The conclusion one may draw from the plots 
is reassuring: while the results for the two bases signifi- 
cantly differ for excessively small sizes [cf. the blue and 
green curves for relative-motion orbitals {Nf, = 2) and 
Slater determinants (two energy shells of Fock-Darwin 
levels), respectively], they tend to substantially coincide 
if the basis size is sufficiently large (e.g. iV^ > 4 or more 
than four energy shells). Besides, a rule of thumb may 
be stated which connects the two bases: The CI diago- 
nalization allows for essentialy the same results using ei- 
ther the Slater-determinant or the relative-motion bases, 
provided Nt equals the number of energy shells of Fock- 
Darwin levels which are filled in to build the Slater deter- 
minants. This is seen in Figs. lllf a) and (b), as the curves 
for the Slater-determinant and relative-motion bases — 
if matched according to the rule of thumb — tend not 
only to overlap as the basis sizes increase, but also to 
present the same dependence from g. The latter is given 
by the separation between two consecutive circles deco- 
rating the curves of Figs. [TlTa) and fTTT b). corresponding 
to the increment of the coupling constant Ag = 1. 

The above rule of thumb may be explained as follows. 
The fact that the number of energy shells of Fock-Darwin 
levels must at least equal Nf, is a necessary condition. In 
fact, the ground-state wave function, which is factorized 
into r- and i?-dependent factors in the relative-motion 
frame, in the Slater-determinant space is built as super- 
position of products of two Fock-Darwin orbitals which 
depend on ri and r2, respectively. Since Laguerre poly- 
nomials up to order Nj, — 1, LArj,_i(r^/2), span the r- 
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FIG. 12: (Colour) CI energy of the two-fermion ground state, 
E, vs. coupling constant, <; > (a) and g < Q (b), for dif- 
ferent numbers of energy shells of Fock-Darwin levels used 
in the CI expansion of the state over a basis of Slater de- 
terminants. The vertical error bars are estimated according 
to the procedure described in the main text. The right ver- 
tical axis shows the two-dimensional scattering length, a2D, 
obtained through Eq. (|10p . The center-of-mass contribution 
has been subtracted from the total energy, therefore the non- 
interacting energy is _E = 1. 



dependent factor in the relative-motion frame, and be- 
cause LAr^_i(r^/2) is a polynomial of order iV;, — 1 in 
r^, then the ri- and r2-dependent Fock-Darwin orbitals 



must include powers of the radii, r^'^^ and r^^^ , respec- 
tively, at least of order f3i+ (}2 = N^ — 1. By considering 
only Fock-Darwin orbitals with zero azimuthal quantum 
number, m = 0, which typically provide the strongest 
contribution to the singlet ground state [s^, the above 
constraint is enforced only if the number of energy shells 
is not smaller than Ni,. With regards to the sufficiency 
of the rule of thumb, we note that, as the number of 
energy shells increases, the neglected high-energy con- 
tributions to the sum which gives ^(r) gradually lose 
their relevance. Therefore, the usage of the two basis 
sets becomes essentially equivalent. Note that this same 
argument holds also for two bosons. 



For the sake of completeness, in Figs. [T2r a) and (b) 
we show the dependence of E (or 020) on g as com- 
puted from a CI diagonalization over a basis of Slater 
determinants. These curves should be directly compared 
with Figs. fTOT a) and (b), obtained in the relative-motion 
framework. Since the plots in Fig. [T^] provide informa- 
tion regarding the estimated error (vertical bars) asso- 
ciated to a given size of the Hilbert subspace (number 
of energy shells of Fock-Darwin levels) , they are particu- 
larly useful for actual CI calculations. The detailed com- 
parison between Fig. [H^a) [Fig. [H^b)] and Fig. 
[Fig.fTUTb)] confirms the rule of thumb: the curves are sig- 
nificantly different for small sizes of the diagonalization 
space [cf. the orange curves in Fig. fT^ a) and Fig. [TUTa)]. 
while they tend to coincide as the size of the space in- 
creases, especially for repulsive interactions. For exam- 
ple, the asymptotic values of the blue, red, and black 
curves in Figs. [T^ a) and fTUf a) differ only for a few parts 
per thousand, and the differences become smaller as the 
basis size increases (from blue to black). 



VIII. CONCLUSION 

The aim of this paper is to provide a ground for CI 
calculations of TV particles interacting via a contact po- 
tential. The configuration interaction method cannot 
be naively applied in two and three dimensions, due 
to the apparent pathology of the (5-function. To cure 
the problem for practical applications in a very simple 
and straightforward manner, we developed a procedure 
to renormalize the strength of the contact potential for 
two particles in a given subspace of single-particle basis 
states. The final outcome is that the CI diagonalization 
over a finite basis set provides physical observables whose 
values do not depend on the energy cut-off of the basis. 

This procedure relies on the comparison of both the en- 
ergy and wave function of two particles obtained by two 
different sets of data. The first set is computed by CI di- 
agonalization with an energy cut-off, while the second one 
by the exact solution of the Schrodinger equation for the 
regularized form of the contact pseudopotential. Both 
ground and excited states are considered in this compar- 
ison. The analysis of CI data for a truncated Hilbert 
space provides a fully consistent physical picture of the 
results as well as a systematic assessment of the error of 
the calculation. 

We believe this work allows for the application of the 
CI technique to the few-body problem with contact in- 
teractions. The interaction strength and model space of 
single-particle orbitals is defined from the investigation 
of the two-body problem and subsequently used for the 
A^-body problem (see Ref. [Tsl ) . It is natural to assume 
that the mapping between coupling constant, scattering 
length, and size of the Hilbert space for two particles, as 
reported in Figs. [TUf a) and (b), as well as in Figs. [TWa) 
and (b), remains valid for N > 2, since the interaction af- 
fects pairs of particles while scattering events concerning 
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three and more particles are rare for sufficiently dilute 
systems. Besides, as N increases, the number of inter- 
acting pairs, ~ iV^, scales with N differently from the 
number of Slater determinants, « (provided the set 
of single-particle orbitals is fixed) We may therefore 
reasonably assume that our renormalization scheme is 
sound for iV > 2, provided that the set of single-particle 
orbitals is sufRciently large. Nevertheless, if the set of 
orbitals is too small, the A^-body Hilbert space could be 
artificially distorted (space filling effect). 

In Rcf. we presented a first application of the 
method by studying the shell structure and the pairing 
of a few fermions in a two-dimensional trap. We expect 
that the present study paves the way to the ab initio 
investigation of few-body physics in cold atom traps. 
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APPENDIX A: SCATTERING AMPLITUDE FOR 
THE TWO-DIMENSIONAL REGULARIZED 
PSEUDOPOTENTIAL 

In this appendix we derive the asymptotic form of the 
two-dimensional continuum states scattered by the regu- 
larized pseudopotential of Eq. ([S]), V^scudo, in free space. 
In the center-of-mass framework, the problem is governed 
by the Hamiltonian 



free 



pscudo: 



(Al) 



where, in comparison with Eq. ([7]), we have suppressed 
the confinement potential in the Hamiltonian. The scat- 
tering states with positive energy, fc^, are solutions ev- 
erywhere of the eigenvalue problem 



i?,X«'(r) = fc"*(r). 
For /cr 3> 1 the states must have the following form: 



(A2) 



vl/+(r)«e*--/(fc,(^) 



Sirkr 



ikr 



(A3) 



where '^^{f) is the outgoing wave of momentum fc, 
f{k,ip) is the scattering amplitude, defined according to 
the conventions of Ref. 3/, and (p is the azimuthal angle. 
In order to determine the unknow function f{k,(p), we 



parallel Ref. |l7| and write the Lippmann-Schwinger equa- 
tion for ^'^(t') in terms of the two-dimensional Green 
function, 

vl/+(r.) = 6^'=- + Jdr'G+ir, r'; k')V^seudoir')^tir'). 

(A4) 

Here the retarded Green function G+(r, r'; k'^) obeys the 
equation 

(V^ + fc2) G+(r, r'; fc^) ^ S{r - r'), (A5) 
whose explicit solution is 



G+{r,r';k')^-^4'\k\r-r'\), 



(A6) 



where Hq^\z) is the zero-order Hankel function of the 
first kind [3^ . By replacing G"*" and Vpseudo in (|A4p with 
their explicit expressions (|A6p and ([6]), respectively, and 
then performing the integral over r', one obtains 



2 log (^aaoA) 
X \E'+(r')-log(ylAr'^-' ^ 

As the factor 



4'\kr) 



dr' 



c = 



*+(r') -log(y^Ar')r' 



dr' 



(AT) 



(A8) 



in Eq. (|A7p does not depend explicitely on r, we find that 
^'^(r) has the form 



Kir) 



tC 



2 log iAa2DA) 



(A9) 



By inserting Eq. (jA9p into (|A8p and then by explicitely 
performing both d/dr' and lim,p/^o+i after some manip- 
ulation we obtain an explicit expression for C: 



(7 = 



log (ylazpA) 

log (fca2D) 



(AlO) 



By inserting the above value of C back into Eq. (jA9p . we 
obtain an explicit formula for 'l'^(r) which holds every- 
where in space, due to the zero range of the pseudopo- 
tential, and that remarkably does not depend on A: 



Kir) = e}^- 



21og(fca2D) 



U^^\kr). 



(All) 



Finally, by using the asymptotic form of the Hankel's 
function in Eq. (jAlip . 



H^\kr) 



,(fer-V4) forfcr>l, (A12) 



we deduce the expression of /(/c, Lp) by direct comparison 
with Eq. 



HKv) 



27r 



log(fca2D) 



(A13) 
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The isotropic scattering amplitude (jA13p coincides with rameter A, the regularized pseudopotential Ipseudo cor- 
the exact low-energy limit of /(fc, (p) for any short-range rectly reproduces the zero-range features of the true po- 
potential [s^l- Therefore, no matter the value of the pa- tential. 
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